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Abstract 

The evolution and propagation of a partially coherent matter wave 
(PCMW) is investigated theoretically by the correlation function method. 
The ABCD matrix formalism previously used for a fully coherent mat- 
ter wave is extended to the PCMW domain. A new ABCD law is 
derived, using a tensor method to describe the evolution of a PCMW. 
An analytical solution of the first-order correlation function is obtained 
that makes the propagation and evolution of a PCMW very simple and 
clear. As an example, the propagation of a PCMW in a gravitational 
field is calculated numerically. 

1 Introduction 

Matter-wave optics is a new and fascinating branch of research benefiting 
from the rapid progress of techniques for laser cooling and trapping of neu- 
tral atoms and molecules. The achievement of Bose-Einstein condensation 
(BEC) in dilute atomic and molecular gases enables the study of matter 
waves from incoherent to coherent domains pQ. In recent years, coherent 
matter waves have been studied extensively, such as the formation of vortices 
in BEC, linear and nonlinear propagation of atomic laser beams (HG2E], and 
other phenomena. The matter wave has unique applications, as in atomic 
interferometry [5J, atomic holography [SI [7], and so on. 

According to existing theory, pure BEC exists only at zero temperature. 
However, zero temperature is not reachable experimentally due to the limits 
of laser cooling [8j. All the BEC and atomic lasers employed in experiments 
have distributions of momenta [2]. Therefore, matter waves, in the strict 
sense, are not completely coherent. All the cold atom gases produced by laser 
cooling at finite temperature are partially coherent matter waves (PCMW), 
whose properties during propagation and evolution are important issues and 
play a crucial role in the applications of matter waves. Nevertheless, they 
not been investigated sufficiently 
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It is well known that a coherent matter wave can be described by a sin- 
gle wave function whose evolution is governed by the Gross Pitaevskii (GP) 
equation [2]. On the other hand, a PCMW can't be described by a single- 
particle wave function, but requires a much more complicated treatment. 
One method is the finite temperature BEC theory, which separates the mat- 
ter wave field operator into two terms; one term describes the condensate 
part, and the other term deals with the noncondensate part. However, the 
equations based on this method are complicated and can't be solved ana- 
lytically. Additional approximations are necessary to solve the equations. 
For example, in the Hartree-Fock-Bogliubov (HFB) approximation [9], the 
condensate wave function satisfies the GP equation, and the noncondensate 
operator is described by an equation that follows from subtracting the GP 
equation from the total field operator's Heisenberg equation. In the Hartree- 
Fock (HF) and Popov approximation [10], the condensate wave function is 
defined by the modified GP equation, whereby the "anomalous" density is 
neglected and the noncondensate part consists of thermally excited atoms as 
described in terms of a semiclassical phase space distribution function. In an 
extended HFB theory, the collisionless noncondensate dynamics is included 
within second-order perturbation schemes, and the expressions for damping 
rates and frequency shifts of low-energy modes are derived [Tl] . 

Another approach to describe a PCMW is the correlation function, which 
provides quantitative information about the coherence and the intensity of 
matter waves. The first-order correlation function characterizes local fluctu- 
ations of the phase of the complex matter wave field amplitude, and is related 
to the contrast achievable in an interference experiment. The second-order 
correlation function is related to fluctuations of the modulus of the complex 
matter wave field amplitude, and expresses the tendency of atoms either to 
cluster or to remain spatially separated [12]. Meystre et al. derived the Van- 
Citter-Zernike theorem of PCMW, which disregards the interactions among 
the particles. They analyzed the propagation of the matter wave from an 
incoherent matter wave source p3], and discussed the detection of the cor- 
relation function [14j . This method is analogous to that used in traditional 
optics [151 HU C£7] . In this paper, we focus on the propagation and evolution 
of PCMW by applying the correlation function method, and disregard the 
interactions among the particles, although it is not necessary. The ABCD 
matrix formalism [181 19] , which previously had only been applicable for the 
wave function of a single particle, is extended to the PCMW. A generalized 
ABCD law to describe the evolution of the PCMW is derived by a tensor 
method. An analytical solution of the first-order correlation function after 
evolution is obtained, which makes the evolution problem of the PCMW 
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very simple and clear. The results are useful in the analysis of the space- 
time coherence, spatial distribution, propagation and imaging properties of 
the PCMW. 

This paper is organized as follows: In Secj2l we introduce briefly the 
ABCD matrix formalism for a single-particle wavepacket. In SecJSl the 
ABCD matrix formalism is extended by a tensor method to describe the 
PCMW. In SecHl the evolution properties of an ultracold atomic sample 
released in a gravitational field are calculated numerically and discussed. A 
summary is given in the Conclusion. 

2 The ABCD matrix formalism 

First, let us briefly describe the ABCD matrix formalism for a matter wave, 
introduced by Ch. J. Borde based on the time-dependent Schrodinger equa- 
tion [IB], which can be used to calculate the evolution of a single atom 
wavepacket \ip(t) >: 

«9*!a> = H\m > ■ (i) 

This equation can be solved in terms of the quantum propagator [18], which 
can be obtained by the use of a shortcut through the classical limit, and 
a well-known result of Van Vleck that makes the connection with quantum 
mechanics. 

The Hamiltonian of a single atom with a mass m in a gravito-inertial 
field in the classical limit is 

H = p T /3p/(2m) - n T L - mg T q 

-mq r 7q/2 + y(t), (2) 

where V(t) is the external field, q T = (x,y,z) is the position vector, p T = 
{PxiPyiPz) is the momentum vector, and L is the angular momentum vector 
of the atom. The gravitational wave and gravitational gradient are repre- 
sented by the tensors (3 and 7. The quantity fl is the angular velocity vector 
of the Earth's rotation, and g is the gravitational acceleration vector. The 
Earth's gravito-inertial field is characterized by the above four parameters: 
(3, 7, Q, g. The superscript T means the transpose. Some judiciously chosen 
transformations and choices lead to the following ABCD formalism for the 
evolution of the variables q and p from the initial time t a to the final time 

t b m- 

( Pb/m ) = ( n^Ui ) + ( C D ) ( pjm ) ' (3) 
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where n _1 is the matrix corresponding to the gravitational wave tensor f3. 
Were we to neglect the gravitational wave effect, then n^ 1 would be a unit 
matrix. The matrix U corresponds to the rotation term —fl T L in Eq.(f2|). 
If there is no angular momentum of the atom, U is also a unit matrix. 
The vector £ is the displacement produced by gravity, which describes the 
classical trajectory of a non-rotating Hamiltonian with the initial conditions 
£(ta) = and £(t a ) = 0. The 3x3 matrices A, B,C,D and the vector £ 
can be determined by solving the Hamilton- Jacobi equation|18j. 

From the above expressions, by neglecting the gravitational wave effect 
and the rotation, the classical action can be written as 

S (q&,£&><la)£a) 

= mi ■ (q& - £) + I " L^t^dt! - f " V(t 2 )dt 2 

J t a Jt a 

+y {(<£ ~ ODB-^qt, - - 2(q£ - ^ T )(B- 1 ) T q a 
+q^B- 1 Aq a ], (4) 

where L\(t) is a partial Lagrangian given by L\(t) = m(|£| 2 + £ T 7£ + 
2g T £)/2. Knowing the classical action allows the determination of the quan- 
tum propagator according to Van Vleck's formula[20j: 

K(q b ,t b ,q a ,t a ) = {—) \detB\~ 1/2 exp(iS bta /h), (5) 

where h is the Plank constant and H = hj (2ir). A complete set of solutions of 
the Schrodinger equation can be derived from K. The lowest order Gaussian 
wavepacket at initial time t a is given by 

i/j(q,t a ) = ip wp (q,t a ;q a ,v a ,x a , Y a ) 

' ^xpS^-qDY.X^q-qJ] 



Vdet X a L 2h 
xexp[^v a - (q-q )], (6) 

which is centered at position q a , has an average velocity v a = p a /m, and 
the complex width parameters X a , Y a . Its evolution 

tp(q, t b ) = J K(q, t b , q', t a )ip wp {q' ', t a ;q a , v a , X a , Y a )dq' 

= exp[ ^%g^ ]^ wp (q, t a ; q b , v b , X b , Y b ) (7) 
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is governed by the ABCD matrix formalism for q and p in Eq.([3|) and for 
X and Y: 

X b = AX a + BY a , (8) 
Y b =CX a + DY a . (9) 

From the ABCD matrix formalism, the problem of the evolution of a 
single particle wavepacket can be solved by calculating the matrices A, B, C 
and D, from which the classical action and the classical equations of motion 
can be obtained. 

As a many-particle system, such as the PCMW, can't be characterized 
by a single-particle wavepacket, the theory must be extended. A convenient 
approach to this problem is to carry out a second quantization, which will 
be done in the following section. 



3 Tensor ABCD law for the PCMW 

In second quantization, atoms can be described by a quantum field operator 
V>(q, t) that satisfies the commutation relations 

$(<l,t),ft(q[,t)]± = tf(q-q'), (10) 
$(q,t),^(q',*)]± = 0, (11) 

where t/m is the conjugate operator of ip. The expression [• • •]_ is a commu- 
tator for bosons, while [• • •]+ is an anticommutator for fermions. The field 
operator ?/>(q, t) is interpreted as annihilating a particle at (q, t), and V^(q, t) 
represents the creation of a particle at (q, t) . The physical properties of a 
many-atom gas can be expressed in terms of correlation functions that are 
expectation values of the field operators, such as the first-order correlation 
function 

r(qi,ti,q 2 ,t 2 ) = (^(Ql, W(q2,t 2 )>, (12) 

and its higher orders: ($'(qi, ti)'0T(q 2j t 2 ) ■■■ ^(q-j) ^^(qij H ere we 
shall restrict ourselves to the first-order correlation function of bosonic par- 
ticles in the form of Eq. (|12p . Its diagonal element I\q, t,q, t) is the atom 
number density at position q and time t. 

For a completely coherent matter wave, the first-order correlation func- 
tion can be factorized into the form [12] 

(^( qi ,il)^(q 2 ,i 2 )) = V*(qi,tiMq2,t2), (13) 
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with 

^(q,t) = VN(q\<Kt)), (14) 

where atoms occupy the same one-particle state \4>(t)), and N is the total 
atom number. The degree of coherence 

, , . s r(qi,ti,q 2 ,t 2 ) nK s 

<?(qi, *i, 02, i 2 ) = , = , == (15) 

vr(qi , ti , qi ,*i ) v r (q 2 , *2 , q2 , *2 ) 

is unity for a perfectly coherent matter wave such as the pure BEC, which 
contains millions of atoms associated with a single-particle wave function. 

The ABCD matrix formalism in SecE] deals with the evolution of a one- 
particle wavepacket and is appropriate to describe a coherent matter wave 
disregarding interactions among atoms [211 \22\ l23j [24] . If the coherence of 
the matter wave is not perfect, that is, if 

0< 5 (qi,ti,q 2 ,t 2 ) < 1, (16) 

it is not possible to factorize the first-order correlation function into the 
form of Eq. (113j) . That means that the matter wave field is only partially 
coherent and some randomness exists. The field can't just be represented 
by the single-particle wave function, and second quantized field theory has 
to be applied. 

We concentrate on the propagation and evolution of non-interacting mat- 
ter waves from a PCMW source. Neglecting the gravitational-wave effects 
and rotational effects, the single-particle Hamiltonian of Eq.([2]) turns out to 
be 

H = p T p/(2m) — mg T q — mq- r 7q/2. (17) 

Changing this single-particle Hamiltonian to the second quantized Hamil- 
tonian, the Heisenberg equation of motion for the field operator ^>(q, t) can 
be found [25J: 

„a *. . / h 2 X7 2 rp m \ 
ih—v(;{q, t) = I — ~ m § q ~ T q 7q J *)■ ( 18 ) 

It has the same form as the Schrodinger equation for the single-particle wave 
function, so solution techniques are similar. The propagation of a matter 
wave field is described by the quantum mechanical propagator 

$(q, *) = J K (<li t, qo, *o)^(qo, t Q )dq . (19) 
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The action Eq.(JH) can be rewritten in tensor form as 
S (q,i,q ,t ) 



™£-(q-£)+ / Li(h)dti 

T 



to 



V(t 2 )dt 2 



to 



m 



qo 

q-£ 



B *A -B" 1 
C-DB *A DB 1 



qo 

q-l 



(20) 



From Eqs.©, (|12p and (|19p . we find that the general propagation formula 
for the first-order correlation function is 



r(r,t) 



1 



detBl 



r(r ,i )exp 



.AS\ 



(21) 



where r and ro are position tensors given by r T = (q^,q^), and Tq = 
(Qoi 5*302)) = (Si — S 2 ) is the action difference, where Sj(J = 1,2) is 
the action from the incident point qoj to the output point qj, and can be 
obtained from Eq. ([20l) by replacing q with qj, qo with qoj. The action 
difference is 



AS = mi ■ (qi - q 2 ) + 



in 



r ° _ 



T 



V 



r ° _ 



where £ T = (£ T ,£ T ), 



V 



B~ 1 A -B 
C-DB 1 A DB 1 



and 



A 
C 



A 
A 

C 

-c 





B 


, B = 


-B 





D 


" 


, D = 







D 



Since AS is a scalar, we get the relations 



(B^A) 



B _1 A, (DB M T 



DB 



C — DB A 



(B) J , 



(22) 

(23) 

(24) 
(25) 



(26) 
(27) 



among the sub-matrices of Eq. (|23p . These will be used in derivations that 
follow. 
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In conventional optics, the cross-spectral density of a partially coherent 
wave is usually described by the Gaussian-Schell-model(GSM)[l5j expression 

W(p 1 ,p 2 ) = ^/S{pT)^S{p~ 2 }g(p 1 ,p 2 ), (28) 

where S(p) represents the spectral density, and g(pi, p 2 ) represents the spec- 
tral degree of coherence, given by 

S(p) = G eM-\p T {<7 2 s r l p\, (29) 

g(pi,p2) = exp[--(pi - p 2 ) T (<r 2 g y 1 (p 1 - p 2 )], (30) 

where 

Pi = (xi, yi), p\ = (x 2 , y 2 ) (31) 

are position vectors of two arbitrary points in the transverse plane. Go is a 
positive quantity. cr s represents the transverse spot width and a g represents 
the transverse coherent width. cr s and cr g are 2x2 matrices with transpose 
symmetry. 

This partial coherence theory in optics can be extended to describe the 
PCMW according to the analogy between the Schrodinger equation (fT8l) of 
a matter wave and the paraxial wave propagation equation of light 

V 2 T A(q) + 2ik^^ = 0, (32) 
oz 

where A(q) is the amplitude of the electric field as expressed by J5(q, t) = 
A(q) exp[i(kz — ut)], and k = 2ir/\ is the modulus of the wave vector. The 
squared transverse Laplacian is mathematically analogous to the atomic 
kinetic energy — (fo 2 /2m)V 2 . The difference between Eq. (|18p and Eq. (j32[) is 
that Eq. (|18p has a time derivative while the paraxial wave equation Eq. (|32p 
has a spatial derivative along the propagation axis z. The GSM partially 
coherent light as described by Eq.( I28p has two spatial dimensions perpen- 
dicular to the propagation axis z, while the GSM of the PCMW may have 
three spatial dimensions. The first-order correlation function of a GSM 
matter wave can be expressed as 

r(qi,q 2 ) = G exp{ - \ [qR^J^qi + qi'W^qa] 

-^(qi-q 2 ) T K)- 1 (qi-q 2 )}, (33) 

where both the atom number density and the coherence degree have Gaus- 
sian distributions. The parameters Go, cr s and cr g are related to the tem- 
perature of the atom gas. 
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Equation (|33j) can be rewritten as 



r (r) - (l - |(lxI) (^M7i r 



(34) 



with r T = (qi,q2)- The tensor 



m; 



2m 



2m' 



(35) 



is a 6 x 6 symmetric matrix, which may be called the complex curvature 
tensor of the PCMW; where cr s and cr g are 3x3 matrices with transpose 
symmetry, given by 





O sxx 


@ sxy 


O sxz 


<?s = 


O syx 


&syy 


&syz 




O szx 


Q szy 


O ' szz 




&gxx 


&gxy 


Vgxz 


a 9 = 


Vgyx 


a gyy 


a gyz 




@gzx 


°gzy 


® gzz 



(36) 



(37) 



The tensor expression can be used to describe a general matter wave system 
including the asymmetric and anisotropic systems, and is very convenient 
for its compactness. The matrix cr s , representing widths of the matter wave 
in 3D space, characterizes the spatial size of the matter wave. The matrix 
(T g describes the coherent length of the matter wave. When a g decreases 
to zero, it corresponds to a completely incoherent matter wave. When a g 
increases to infinity, the matter wave can be treated as a completely coherent 
matter wave. 

Substituting Eq.J22]) and Eq.flM]) into Eq.{2U, w e get 



r(r,t) 



Go 



m\ 3 

ih) IdetBl 



cxp 



vm 
2k 



Q + 2£ • ( Ql - q 2 ; 



dr , 



(38) 



where 

Q 



(r - 0' (C + DM7 1 ) (A + BMry^r - £) 



+ 



(B *A + Mr 1 ) 1 / 2 ^ - (B^A + Mr 1 )- 1 ^- 1 ^ - £) .(39) 
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In the derivation of Eq. (|39p . the transpose-symmetric property of Mj 1 and 
Eg. ([26]) have been used. Integrating Eq. ([38|) . we have 



r(r) 



G [det(A + BMr 1 )]- 1 / 2 exp 
im 



im£ ■ (q 2 - qi)/h 



x exp 



2h 



(v-tfMjHv-i) 



(40) 



where 1 and 1 denote the complex curvature tensors of the PCMW 
at the initial and the final time, respectively. They satisfy the condition 



M 



-i 
/ 



(C + DMr 1 )(A + BM 



(41) 



In the derivation of Eq. (|40p . the integral formula exp[— ax 2 \dx = \Jtt /a 
has been used. Equation (|41|) may be called the tensor ABCD law for the 
PCMW. 



4 Evolution of the PCMW in a gravitational field 

To illustrate the usage of the formulas derived above, we are going to calcu- 
late the atom density distribution of the PCMW in a transverse plane in a 
gravitational field, which is the diagonal element of the first-order correla- 
tion function T(q, q). The object we consider is a matter wave of cold Rb 87 
atoms. It propagates along the z direction, and the gravitational accelera- 
tion g with a linear gradient 7 goes in the same direction. As is the case 
in Ref . |22j . the cold atom gas is prepared in a magnetic trap with trapping 
frequencies of u x = uj z = 2tt x 330 Hz and uj y = 2tt x 8 Hz. The initial 
temperature of the Rb s7 atom gas is assumed to be the recoil temperature 
by laser cooling [8]: T rc = 362 nk. 

The propagation matrices A, B, C and D are 



(42) 








^ 




(i 


1 





!■»-( 







A z J 













(: 





°) 


...( 







c z ) 





t 


-t 













t-t 



























1 


l \ 










D z J 





(43) 



where A z = D z = cosh[^/7(t-t )], B z = sinh[^/7(t-t )]/V7 1 and C z = -yB z 
The gravitational displacement vector is £ T = (0,0,^), and £ z = (g/j)[l - 
cosh(V7(* - *o))] M- 
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Assuming the ultracold atoms are initially in a magnetic quadratic trap, 
the widths of the atom gas in the x,y,z directions are [2] 

/ \ 1/2 

= (—2) , (44) 

where Uj (j = x, y, z) is the harmonic oscillator frequency in each dimension. 
Substituting the initial temperature T rc and the initial trapping frequencies 
Uj into the expression of a s j, the widths of the matter wave source can be 
obtained: a sxx = 2.84/xm, a syy = 117.25/mi. As we have assumed that the 
ultracold atom gas propagates continuously along the z direction, the width 
along the z axis is infinite: a szz = 00. Then we can write out the width 
tensor of the matter wave source as 

/ \ 

{al)- 1 = a- y 2 mm" 2 . (45) 

V ooy 

We assume the coherent length l c to be IOOAt, which is between the thermal 
de Borgile wavelength At = \/2irh 2 / (uiKbT) and the infinity. The coherent 
length matrix is 

/ i- 2 \ 
(a 2 )' 1 = I' 2 mm" 2 . (46) 
\ / 

Substituting the width matrix (er 2 ) -1 and the coherent length matrix (cr 2 ) -1 
into Eq. (|35p . we can calculate the final complex curvature tensor MJ 1 from 
Eq. (|4*T|) and get the atom density distributions from Eq. ([4"0~|) . The expression 
of A,B, C and D in Eq. (|42p shows that the expansion of the matter wave 
is dependent on time, and the transverse density distributions at different 
temporal points reveal the evolution of the matter wave as it propagates. 

Figure Q] shows the evolution of the PCMW, with the parameters given 
in the caption. The initial atom density profile, which is shown in Fig[T](a), 
is elliptical due to the initial inhomogeneous confinement of the cold atom 
gas on the x-y plane; that is, the confinement in the x-direction is stronger, 
which causes a larger speed of expansion of the matter wave in this direction 
when the confinement is removed. During evolution without confinement, 
the atom density profile expands gradually ( FigHJb)), and changes to a 
circle (FigJTJc)). Finally the atom density profile becomes an ellipse again 
with the major axis perpendicular to that of the initial ellipse, as shown in 
FigETJd). This is because the degree of expansion in the x-direction exceeds 
that in the y-direction. This property is similar to that of light beams. 
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Figure 1: The transverse relative atom density distribution for an anisotropic 
GSM PCMW at different propagating times t. (a) t = Os; (b) t = 0.5s; (c) 
t = Is; (d) t = 2s. 
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5 Conclusion 



In this paper, a theoretical description of a partially coherent matter wave 
(PCMW) based on the analogy between a matter wave and an optical wave 
is presented. The propagation and evolution formula of the Gaussian-Schell- 
Model PCMW is provided, and a generalized tensor ABCD law is derived. 
As an example, we analyzed the atom density profile in a transverse plane 
of a Rb 87 cold atom beam. The results show that the tensor ABCD law is a 
very convenient method for treating the propagation and transformation of 
partially coherent cold atom beams. Most importantly, the previous results 
for partially coherent light can be converted into a PCMW. Our method can 
be applied to analyze the propagation and evolution properties of ultracold 
atom beams in various potentials, such as an optical potential and a mag- 
netic potential. This method can also be extended to pulsed matter wave 
propagation, and to include the interactions among atoms and the internal 
state of atoms with some modifications. These will be the topics of further 
investigations. 
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